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OF LIQUIDS IN TOROIDAL TANKS 

By John Locke McCarty, H. Wayne Leonard, and 
William C. Walton, Jr. 



SUMMARY 


Several toroidal configurations applicable to missile and space- 
vehicle liquid storage systems were oscillated to study the natural 
frequencies of the antisymmetric modes of contained liquids over a 
range of liquid depths and tank sizes. Natural frequencies were obtained 
for tank oscillations parallel to the free surface of both vertical and 
horizontal tank orientations. 

The data are presented in terms of dimensionless parameters which 
are obtained by relating experimentally determined natural liquid fre- 
quencies to analytical expressions developed through consideration of 
the physics of the problem and from existing solutions for liquids in 
tanks having similar boundaries at the liquid surface. The experimental 
results obtained for the toroids indicate that these parameters are 
applicable to the prediction of the natural frequencies of fluids in 
toroids of general geometry and size. 


INTRODUCTION 


In missiles or space vehicles employing liquid- fueled propulsion 
systems or large volumes of liquids for life support, the responses of 
the system to motions of the contained liquids may greatly affect the 
dynamic stability of the entire vehicle. The magnitudes of these 
responses can increase greatly if the natural liquid frequencies are near 
frequencies of external periodic forces which may be induced by control 
impulses or by periodic structural deformations of the vehicle. It is 
imperative, therefore, that methods for accurately predicting the natural 
frequencies of these liquids in propellant tank configurations of interest 
be available before the tank is incorporated into any vehicle design. 

The natural frequencies of the antisymmetric modes are of particular 
interest since the liquid motions in these modes involve lateral shifts 
in the liquid center of gravity. Some results of experimental and 
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analytical studies of the natural frequencies and mode shapes of liquids 
in spheres and right circular cylinders are given in reference 1. Recent 
studies, however, indicate the desirability of employing toroidal storage 
tanks for propellants or liquids for life- support systems in some future 
vehicles, but no information on the natural frequencies of liquids in 
such tanks is presented in the available literature. 

The purpose of this paper is to report the results of an experi- 
mental investigation of the natural frequencies and mode shapes of the 
antisymmetric modes of liquids contained in toroidal tanks. The fre- 
quencies and mode shapes were obtained for tanks having various major 
and minor radii, liquid depths ranging from empty to full, and different 
orientations of the tanks with respect to the direction of the applied 
oscillation. The natural frequencies are presented as dimensionless 
parameters to permit application of the results to toroidal tanks of 
practical interest. 

SYMBOLS 


g acceleration due to gravity 

h liquid depth 

h c liquid depth in an annular right circular cylinder having the 

same liquid- surface geometry and yielding liquid volume 
equivalent to that contained in a horizontal toroid at liquid 
depth h (fig. 7) 

Jp' first derivative of the Bessel function of the first order and 

first kind 

n mode of liquid oscillation 

r minor radius of toroid (fig. l) 

R major radius of toroid (fig. l) 

r^, r Q inner and outer radii of liquid surface for liquids in hori- 
zontal toroids (fig. 7) 

Yj_' first derivative of the Bessel function of the first order 

and second kind 

e n nth zero of Jj_' 
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An 


'n 


■'n 


0 

*n 


angle measured from the vertical to the radius of length R 
that terminates at the liquid surface of vertical tanks 
(fig. 8) 

frequency parameter for nth mode of a liquid in a horizontal 



frequency parameter for the nth tr ansve rse mode of a liquid 
in a vertical toroidal tank, co^r/g 

angle between the vertical and the radius to the point of 
intersection of the liquid surface with the peripheral 
circle of vertical tanks (fig. 8) 


frequency parameter for the nth longitudinal mode of a liquid 
in regions of a vertical toroid (fig. l); for regions A and C, 

* n * % l R * r ( n > i) 

and for region B, 



experimental natural circular frequency of oscillation of 
liquids in nth mode 

analytical natural circular frequency of oscillation of liquids 
in nth mode 


4 


APPARATUS 


Description of Models 

The models of the toroidal tanks studied in these tests consisted 
of five configurations, each of which was oriented in three different 
ways with respect to the direction of oscillation. These orientations 
are shown in the sketch in figure 1, which also includes the tabulated 
dimensions of the models. Orientations of the tanks for the various 
liquid modes were as follows: 


Mode 

Plane of major 

Direction of oscillation with 

torus radius R 

respect to plane of major radius 

Horizontal 

Horizontal 

Parallel 

Longitudinal 

Vertical 

Parallel 

Transverse 

Vertical 

Perpendi cular 


In all modes the direction of oscillation was in a plane parallel to the 
liquid surface. 

All models were constructed of clear Plexiglas to permit visual 
observation of the liquid motion. In all cases water was used as the 
liquid. 


Mechanical Shaker 

The models were mounted on a support platform which was suspended 
in pendulum fashion from overhead beams. Oscillations of the models 
were induced by means of a mechanical shaker which was directly con- 
nected to the support platform as shown in figure 2. The mechanical 
shaker, described fully in reference 2, consists of a slider-crank 
mechanism driven by a variable-speed motor and designed so as to provide 
a means for conveniently varying the frequency and amplitude of the 
reciprocating motion applied to the platform. A tachometer, also shown 
in figure 2, was attached to the drive shaft of the motor and provided 
a means for directly obtaining the excitation frequency. The shaker 
had an additional design feature which permitted a rapid shutdown of 
the driving motion so that the liquid modes excited at a given frequency 
could be studied during the decay of the fluid motions. 
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TEST PROCEDURE 


The testing technique involved inducing translatory oscillations of 
the models over a range of frequencies to obtain the natural frequencies 
of the contained liquid. The procedure was repeated over the full range 
of liquid depths . In measuring the lower mode natural frequencies for 
all mpdels except toroid 5 > the mode in question was induced by means of 
the mechanical shaker and, upon full development of the wave form, the 
platform motion was stopped and the frequencies were obtained by visually 
timing the low-amplitude oscillations of the liquid during the decay of 
the wave form. In the cases of the higher modes, for toroids 1 to 4, 
the excitation amplitudes were maintained at low levels and the natural 
frequencies were taken as those frequencies yielding maximum liquid 
response. In these cases, the frequencies were read directly from the 
tachometer. Inasmuch as toroid 5 was too large for installation on the 
shaker platform, the testing technique for this toroid was modified. The 
toroid was mounted on rollers and manually excited to induce the desired 
liquid mode shape. All liquid modes of this toroid were obtained by 
visually timing the low-amplitude liquid oscillations during the decay 
of the fully developed wave form. Representative mode shapes of the 
liquid in half-filled toroids mounted horizontally and vertically are 
shown in figures 3 to 5- Data were taken for all modes visually detected 
with sufficient clarity for their definition. 


DATA REDUCTION 


A sample of the test results is given in figure 6, which presents 
some of the experimental data taken on toroid 3- This figure shows the 
variation of the measured natural frequencies of the first two liquid 
modes in cycles per second with fullness ratio for the three tank orien- 
tations. In order to apply the results of the experimental tests to 
toroids of variable geometry, an extension of these results to permit 
the prediction of natural liquid frequencies in toroids of different 
dimensions is desired. The ideal method, one using dimensionless param- 
eters which express the ratio of the experimentally determined natural 
liquid frequencies to simple, exact, closed-form solutions for the 
natural liquid frequencies, is not possible at this time since no such 
exact solutions are known to exist. It is possible, however, that cer- 
tain alternate expressions may be derived which, when similarly applied, 
will yield frequency parameters independent of tank dimensions. The 
existence of such expressions is suggested by recognition of the fact 
that at various liquid depths in a toroid, the liquid has physical 
boundaries relatively similar to those of liquids in containers for 
which information concerning natural fluid frequencies is available 
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( e.g ., cylinders, spheres, and U-tubes). For some such containers, 
either simple, exact, closed-form solutions exist for the natural liquid 
frequencies, or a combination of pertinent variables has been proven 
experimentally to yield frequency parameters that are essentially inde- 
pendent of tank dimensions. The derivation of the toroidal frequency 
parameters based upon these expressions is presented in the following 
sections for each mode of oscillation. It should be noted that the fre- 
quency parameters used in the reduction of the experimental data are not 
the only parameters which will nondimensionalize the data. Several 
parameters for each mode were investigated, and those yielding the best 
results were selected for use. 


Horizontal Modes 

The frequency parameter for the horizontal modes was selected as 
the ratio of the experimentally determined liquid frequencies to the 
natural frequencies of liquids contained in an upright annular circular 
cylinder having dimensions r^ and r Q (the inner and outer radii of 
the liquid surface; see sketch in fig. 7) and a liquid depth h c neces- 
sary to produce a liquid volume equal to the volume of liquid contained 
in the toroid. The exact expression for the liquid frequencies 0 ^ in 
such an annular cylinder, given by reference 5, is 



where n is the mode of liquid oscillation, g is acceleration due to 
gravity and v n is the nth root of the equation 



The first four roots of this equation are plotted in figure 7 as a func- 
r • 

tion of — . The resulting parameter for the horizontal modes, denoted 
r o 
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Vertical Modes 

The distinct differences in the geometry of the liquid boundaries 
for various depths in vertical toroids suggest the division of the tanks 
into three separate regions to permit a comparison of the experimental 
data for toroidal tanks of various geometry and size on the same non- 
dimensional fullness basis. These regions are designated A, B, and C 
and are indicated on figure 1. 

Transverse . - In vertical toroidal tanks undergoing transverse oscil- 
lations,' the liquids in regions A (0 < h < 2r) and C (2R < h < 2(R + r)) 
have boundaries which are somewhat similar to those of liquids in spheres 
or horizontal cylinders undergoing transverse oscillations. It has been 
theoretically determined (ref. k) and experimentally verified (ref. l) 
that for a sphere or a horizontal cylinder undergoing transverse oscil- 
lations a frequency parameter cr n of the form 

°„ - “»/! (4) 


will insure satisfactory nondimensionalization of the natural frequencies 
of liquids contained in such tanks. Because of the similarity in the 
liquid boundary conditions, this parameter was selected for use in 
regions A and C. 


In region JB (2r < h < 2R) the liquid boundary is somewhat similar 
to that of liquids in an upright circular cylinder. The exact expres- 
sion for the natural frequencies of liquids in an upright circular 
cylinder of radius a, given in reference 5, is 


^n ~ 



(5) 


where h is the depth of liquid in the cylinder and e n is the nth 
zero of the first derivative of the Bessel function of the first order 
and first kind. In the treatment of the liquid in region B as though 
it were contained in a cylinder of radius r, the reasonable assumption 
is made that the value of the equivalent liquid depth is sufficiently 

large to insure that tanh^ e n ^ approaches unity ^tanh^ -> 1.0 

at ^ = 1.6; tanh^i €2^ -» 1.0 at ^ = 0.6^. With this assumption, 
equation (5) reduces to 


fin - 



(£\ 
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It may be seen from this expression that the prime variables for 
nondimensionalization are g and r - the same variables that con- 
stitute the frequency parameter in regions A and C (eq. (4)). The 
frequency parameter for the upright transverse modes , therefore, is 
made to be consistent throughout the depth range; that is. 


°n = “n 



(7) 


Longitudinal . - The first longitudinal mode of the liquid in 
regions A and C is nondimens ionalized by taking the ratio of the experi- 
mental natural frequency to the natural frequency of a simple pendulum 
whose length is equal to the distance from the center of the toroid to 
the center of mass of a mass distributed uniformly along an arc of the 
peripheral circle subtended by the liquid surface. (See figs. 8(a) and 
8(c).) The expression for the natural frequency of such a pendulum is 


= 


M 


/„ N sin 0 

(R + r) -ir 


( 8 ) 


where 0 is the angle between the vertical and the radius to the point 
of intersection of the fluid surface with the peripheral circle. The 
resulting parameter for the first mode is 


CD, 

+1 , _ = ^ 


R + r sin 0 

1 M d 


(9) 


Nondimensionalization of the higher modes in these two toroidal regions 
was accomplished through the consideration of the relative similarity of 
the liquid boundaries to those of liquids contained in spheres of radius 
R + r. The resulting frequency parameter i|r n is then 


t n = 


\ 


R + r 


(n > 1) 


(10) 


First-mode liquid motion in region B is much like that in the 
circular-arc tube of small cross section illustrated in figure 8(b), 


'l'ne natural irequency oi 


bUC J L 


ovioVi a -hn’hfa is 


0 = I [EeI 

1 \M 9 


(id 
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where R is the radius to the center line of the tube and 0 is the 
angle between the vertical and the radius R which terminates in the 
liquid surface. The dimensionless frequency parameter \|r^ selected 
for the first mode is the ratio of the experimental natural liquid fre- 
quency to the frequency expression given by equation (ll), or 


o>, 

♦i , , ^ 



( 12 ) 


Nondimensionalization of the higher modes is accomplished, as with the 
transverse modes in this region, through the consideration of the rela- 
tive similarity of the liquid boundaries to the boundaries of liquids 
contained in an upright circular cylinder. The order of the analogous 
cylindrical frequency is decreased by 1 since the manometer- type mode 
not present in the upright cylinder was considered as the fundamental 
mode for this region. Again, if it is assumed that the equivalent liquid 

depth is sufficiently large to insure that tanh^ e n j approaches unity 
in the natural-frequency expression (eq. (6)), the frequency parameter is 


= (n>1) <13) 

The liquid frequency parameters for the various toroidal tank 
orientations and the sources for the corresponding analytical expres- 
sions are summarized in table I. 


DATA PRESENTATION AND DISCUSSION OF RESULTS 


The frequency parameters are plotted as a function of fullness 
ratio h/2r for the horizontal toroidal tanks and as a function of 

region fullness ratio ^region A, h/2r; region B, — ~T> re gi° n C, 

ll — p~R \ 

' for the vertical toroidal tanks. Data for each mode 


2r 


of a given 
The data 


configuration are presented in this form in figures 9 "to 11. 
obtained for the higher modes of the small models were limited because 
the mode shapes could not be clearly defined in some cases. 


Modes of Horizontal Toroids 

The experimental, data for the first four modes of liquids in five 
horizontal toroidal tanks are presented in figure 9 in terms of the 
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frequency parameter A n . The values of this frequency parameter are 
plotted as a function of the fullness ratio. The figure indicates that 
at a given fullness ratio, the frequency parameter for a given mode is 
the same for all toroids examined. It appears, therefore, that the fre- 
quency parameter A n is independent of the tank dimensions (both actual 
size and geometric ratio R/r) and applicable for predicting the natural 
liquid frequencies in any horizontal toroidal tank. 

It is of interest that, except at the near-full and near-empty con- 
ditions, the values of A n are in the neighborhood of unity, which 

indicates that the analytical expression for the frequency (eq. (l)) is 
a good approximation to the natural frequency of the liquid in horizontal 
toroids. Both this expression and the frequency parameter derived there- 
from may be readily utilized to obtain the liquid frequencies by using 
the values of v n presented in figure 7 for the first four modes. 


Transverse Modes of Vertical Toroids 

The results obtained from the transverse oscillations of the vertical 
toroids are presented in figure 10. The frequency parameter a n is 
shown as a function of fullness ratio of the three toroidal regions for 
the first two liquid modes. As in the case of the liquid modes in hori- 
zontal toroids, the nondimens ional data obtained for a given transverse 
mode of the vertical toroids may be represented by one curve. It appears, 
therefore, that the resulting curve is applicable for the prediction of 
the natural liquid frequencies in vertical toroids of any size and 
geometry undergoing transverse oscillations. 


Longitudinal Modes of Vertical Toroids 

Liquid frequency-depth relationships for the longitudinal modes of 
vertical toroids are presented in figure 11 in terms of the frequency 
parameter i|r n . The values of this parameter for each of the three regions 
are plotted separately as a function of the region fullness ratio. The 
results for region A are presented in figure 11(a) and indicate that for 
a given fullness ratio, the frequency parameter for the first mode tq 
is the same for all toroids examined. The figure also shows that the 
values of i|r-^ are near unity throughout the region, which indicates 
that the expression given by equation (8) essentially predicts the first 
longitudinal liquid frequency in region A of vertical tanks. 

The frequency parameter developed for the higher liquid modes in 
region A, based upon an expression involving the variables necessary to 
nondimensionalize satisfactorily the transverse modes of spheres and 
horizontal cylinders, is also presented as a function of region fullness 
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ratio in figure ll(a) . The data show that this parameter, (where 

n > l), while dimensionless, is dependent on the ratio R/r, hut inde- 
pendent of tank size. The effect of R/r on the second mode, however, 
appears to he negligible. The faired data of the higher modes show that 
the value of the frequency parameter at a given liquid depth increases 
with this ratio R/r. 

The results of the frequency data obtained in region B are shown in 
figure ll(b) and indicate that the frequencies for the three modes 
examined in this region are readily nondimens ionalized by the parameters 
selected. The natural-frequency parameter i|r n (n = 1, 2, 3 ) has the 
same value for all toroids examined at a given region fullness ratio and 
a given mode. Thus it may be concluded that the parameters are inde- 
pendent of the tank dimensions and applicable in the prediction of the 
longitudinal liquid frequencies in the midregion of vertical toroids of 
practical interest. Furthermore, since the parameters are near unity 
except at the near-full condition, the analytical expressions for the 
frequencies may be used to obtain approximate natural frequencies through- 
out most of the depth range of region B. 

The variation of the frequency parameters \|r n with fullness ratio 
for region C is presented in figure ll(c). As was the case for region A, 
the frequency parameters selected appear to be sufficient to reduce the 
data to a family of curves indicating the dependency of the natural fre- 
quencies on R/r. Thus far, no effective parameter has been found which 
is useful for combining the effects of variations of both tank size and 
geometry in this region. The data for a given geometric ratio may be 
extended directly to geometrically similar tanks. 


CONCLUSIONS 


Frequency data have been obtained and parameters synthesized for 
liquids in toroidal tanks of different sizes and geometric ratios R/r 
(where R is the major radius and r the minor radius of the toroid). 
These parameters are developed for tanks oriented in three ways with 
respect to the direction of oscillation and have been found to be inde- 
pendent of tank size and in most cases independent of the geometric 
ratio R/r. These parameters are believed to be applicable in the pre- 
diction of the natural frequencies of fluids in toroids of practical 
interest in missile and space-vehicle liquid-storage systems. 

The natural frequencies of liquids in horizontal toroids in general 
may be accurately predicted by the analytical expression for the modal 
frequencies of liquids in annular cylinders having the radial dimensions 
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of the toroidal liquid surface and a depth yielding a volume equivalent 
to the volume of liquid in the toroid. 

The transverse modes of liquids in the upper and lower regions of 
vertical toroids behave in a manner similar to modes of liquids in 
spheres, whereas the modes in the midregions behave in a manner analogous 
to modes of liquids in deep upright circular cylinders . 

For the longitudinal modes of liquid in vertical toroids no method 
has been found to achieve complete nondimensionalization with respect 
to both size and geometry, but as expected, it was found in all cases 
that natural- frequency data for geometrically similar toroids can be 
rendered dimensionless and thus the effects of tank size can be isolated. 
The data presented indicate that the degree of dependence of the natural 
frequencies on R/r is a function of both liquid depth and the order of 
the liquid mode. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., July 15, i960. 
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TABLE I.- SUMMARY OF NONDIMENS IONAL FREQUENCY PARAMETERS 
FOR LIQUIDS IN TOROIDAL TANKS 


Orientation 

Region 

Mode 

Parameter 

Source 

Horizontal 

All 

All 

■ 

Axmular circular 
cylinder 

Vertical 

transverse 

B 

All 


Sphere of radius r 

B 

All 


Circular cylinder 
of radius r 

Vertical 

longitudinal 

1 

First 


Simple pendulum 

All n > 1 

(R + r 

Sphere of radius 
R + r 

B 

First 


Circular-arc tube 

All n > 1 


Circular cylinder of 
radius r 


a For convenience, values of are listed below: 

n 

e± = 1.811 
e 2 = 5-331 
= 8.536 
= 11.706 





































Longitudinal excitation 


Transverse excitation 



Toroid 

R,in. 

r,in. 

1 

4.0 

4.0 

2 

7.9 

3.9 

3 

11.9 

5.9 

4 

1 1.9 

3.9 

5 

16.0 

8.0 


Figure 1.- Sketch showing orientation of toroidal tank and dimensions of 

test configurations. 
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L-6o-i8i6.i 

Figure 2.- Test apparatus for mechanically exciting toroidal tanks. 
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(a) First mode. 



(b) Second mode. 



( c ) Third mode . 


L_ 60-4306 


Figure 3*- First three modes of liquid in a horizontal toroid. 


(a) First mode. 



(b) Second mode. L-60-4307 

Figure 4.- First two transverse modes of liquid in a vertical toroid 




(a) First mode. 



(b) Second mode. 



( c ) Third mode . 
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Figure 5-- First three longitudinal modes of liquid in a vertical toroid. 






toroidal tank. (Toroid number 3.) 




(a) Region A. 



( b ) Region B . 



(c) Region C. 


Figure 8.- Sketch defining symbols used in treatment of longitudinal 
liquid modes in vertical toroidal tanks. 
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Fullness ratio 
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Figure 9*- Variation of liquid frequency parameter 


A n « /- 

1/ S Vn tanh 


h„ 

*o Vn i 


with depth for hori 


zontal toroidal tanks . 



Regional fullness ratio 


Figure 10.- Variation of liquid frequency parameter cr n = 

depth for transverse modes of vertical toroidal tanks 


















fullness ratio Regional fullness ratio 



Figure 11.- Concluded. 









